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Abstract 

It has been recently suggested by Dvali and Vikman that the superluminal neutrino 
phenomenology of the OPERA experiment may be due to an environmental feature of the 
Earth, naturally yielding a long-range fifth force of gravitational origin whose coupling with 
the neutrino is set by the scale M*, in units of reduced Planck mass. Its characteristic length 
X should not be smaller than one Earth's radius R e , while its upper bound is expected to 
be slightly smaller than the Earth-Moon distance (60 R e ). We analytically work out some 
orbital effects of a Yukawa-type fifth force for a test particle moving in the modified field of 
a central body. Our results arc quite general since they are not restricted to any particular 
size of X] moreover, they are valid for an arbitrary orbital configuration of the particle, 
i.e. for any value of its eccentricity e. We find that the dimensionlcss strength coupling 
parameter a is constrained to |oc| < 1 x 10~ 10 — 4 x 10~ 9 for 1 R c < X < 10 R e by the laser 
data of the Earth's artificial satellite LAGEOS II, corresponding to M*>4 x 10 9 - 1.6 x 10 10 . 
The Moon perigee allows to obtain |oc| < 3 x 10 -11 for the Earth-Moon pair in the range 
15 R e < X < 60 R e , which translates as M*>3 x 10 10 — 4.5 x 10 10 . Our results are neither 
necessarily limited to the superluminal OPERA scenario nor to the Dvali- Vikman model, 
in which it is M» ~ 10~ 6 at X ~ 1 R e , in contrast with our bounds: they generally extend 
to any theoretical scenario implying a fifth-force of Yukawa-type. 

PACS: 04.50.Kd; 96.30.-t; 95.10.Eg; 14.60.Lm; 14.60.St; 13.15.+g 

1 Introduction 

The measured arrival times 

6t = 57.8 ± 7.8 Ut Tlg| syst ns ~ 58 ± 13 ns (1) 

of the muon neutrinos in the OPERA experiment pQ after a travel along a known baseline 
distance d = 730 km from the CERN CNGS beam to the Gran Sasso LaboratorjQ (LNGS) 

1 See also the material pertaining to the conference held at CERN on September 2011 here: 
http:/ /indico.cern.ch/conferenceDisplay.py?confId=155620. 
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have been interpreted in terms of a property pertaining the motion of the neutrinos themselves. 
Indeed, it was suggested that their speed v„ may have overcome the speed of light in vacuum. 
Stated differently, by defining the dimensionless quantity 

i = P) 

C 

it would be [1] 

Cmcas > (3) 

at a statistical significant level. The result [1] 

£ = (2.37 ± 0.32| stat +^| syst ) x 10- 5 ~ (2.4 ± 0.5) x lO" 5 (4) 
was interpreted by Adam et al. [1] by assuming 

C = <' + A %. ( 5 ) 

with 

A V =c(Ws + l)-< } >0, (6) 

where v^J denotes the expected value of the neutrinos' speed: vffl < c since the neutrinos have 
a tiny non-zero mass. It turned out [1] that the measured effect, within the accuracy of the 
measurement, does not depend on the energy of the muon neutrinos in the domain explored by 
OPERA (some tens of GeV). Such an interpretation of the neutrinic phenomenology observed 
in OPERA is unavoidably bound to raise many questional [21 El II El [6]. However, 0E] pointed 
out that the statistically significative positive result of eq. (jlj) might be due to a systematic 
error of instrumental origin, so that £ would now be statistically compatible with zero within the 
error bar. Recently, a time-of-flight measurement of neutrinos from CERN to LNGS performed 
by the ICARUS collaboration [9j did not confirm the OPERA result by reporting 

St = 0.3 ± 4.0| stat ± 9.0| syst ns ~ 0.3 ± 10 ns, (7) 

which is statistically compatible with zero. Note that, as far as the uncertainty in £ is concerned, 
eq. ([7]) naively corresponds to about the same level of eq. d3J, i.e. about 0.5 x 10 -5 ; indeed, the 
accuracies in determining St are about the same in both eq. ([I]) and eq. 0. For earlier studies 
on superluminal motions of neutrinos in various frameworks, see, e.g., |10 1 111 ! [T2 l 113 1 [T l 115 1 116|. 

Dvali and Vikman [17] wondered if the OPERA superluminal phenomenology could be an 
environmental effect characteristic of the local neighborhood of the Earth, without the need 
of violation of the Poincare invariance at a fundamental level. Such a scenario, at an effective 
field theory level, yields naturally an inevitable appearance of a testable long-range fifth force 
of gravitational type. Following a clarification by Dvali [18], it is 

M* ~ (8) 



In it, M P = = 2.43 x 10 18 GeV is the energy equivalent of the reduced Planck's mass, 

where h is the reduced Planck constant, c is the speed of light in vacuum and G is the Newtonian 



2 In view of the expected forthcoming large amount of papers dealing with [T] , we will not try to formally cite 
them here because such a list would likely become out-of-date very quickly. See the electronic databases like 
ArXiv, NASA/ADS, SPIRES/HEP. 
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constant of gravitation, and the scale M* sets the strength of the coupling of the putative new 
massive spin-2 degree of freedom to the neutrino [17] . For other investigations involving various 
aspects of gravitation, astrophysics and cosmology, see, e.g., P32 EU1 EH E21 E31 [Ml ESI ESI EH 
ESI ESI [301 EI] • According to Dvali [18] , the range length X should not be shorter then the 
terrestrial radius R e : see also [32], in which Earth-size extra dimensions were studies. Data 
from solar neutrinos, not yet analyzed in this respect, would allow to obtain an upper bound 
on X which would likely be shorter than the Earth-Sun distance [18]. Thus, Dvali [18] argues 
that X should be something less than the Earth-Moon distance, and larger than Earth's radius. 

In this paper, we will explore such an intriguing possibility by analytically working out 
some orbital effects of a gravitational long-range fifth force of Yukawa-type. We will perform a 
first-order perturbative calculation without making any a-priori assumptions either on the size 
of X or on the trajectory's configuration of the test particle orbiting the central body acting 
as source of the putative exotic effect. We will also put constraints on the strength parameter 
a for various values of X within the ranges envisaged by Dvali and Vikman [T7] in view of the 
latest results from the orbital determination of some natural and artificial bodies around the 
Earth. Thus, it will be possible to infer lower bounds on the coupling of neutrino to the new 
hypothetical force, which is given by Mp/M* [T7]. The upper bound on Mp/M* comes from 
astrophysical and cosmological observations like star cooling and Big-Bang Nucleosynthesis 
(BBN) [32] : a^l Dvali and Vikman have shown [T7], the requirement that new particle is not 
affecting star-cooling and BBN yields M*/Mp > 4 x 10~ 12 — 10~ n or so, corresponding to 
M* > 10 7 - 10 8 GeV in equivalent energy units. 

The plan of the paper is as follows. In Section EJ we analytically work out some long- 
term orbital effects due to a Yukawa-like modification of the Newtonian inverse-square law. 
In Section 13.11 and Section 13.21 we phenomenologically constrain a with the laser data of the 
Earth's artificial satellite LAGEOS II and the Moon, respectively. We also infer corresponding 
bounds on M*/Mp. Section H] summarizes our findings. 



2 Analytical calculation of some orbital effects induced by a 
Yukawa-like fifth force 



The Yukawa-type correction to the usual Newtonian gravitational potential — ~^/ r i where 
H = GM is the gravitational parameter of the central body of mass M which acts as source of 
the supposedly modified gravitational field, is 



AUy 



■ cxp 



(9) 



where is the gravitational parameter evaluated at distances r much larger than the scale 
length X. The total acceleration resulting from 



Ui 



tot 



C/ N + AUy 



r 



1 + a exp 



is, thus, 



.4 



tot 



/J'OO 



1 + OC 1 + 



X 



exp 



(10) 



(11) 



It should be noticed that Dvali and Vikman [T7] leave room, in principle, for a composition- 
dependent fifth-force, so that a may not be the same for different bodies. From eq. ([lip the 



Dvali and Vikman [17] used the reduced Planck units, in which Mp = 1, and c = 1. 
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following considerations can be traced about the relation between /i^ and the values /i mea s 
of the gravitational parameter actually measured in, e.g., ranging experiments to terrestrial 
artificial and natural satellites, interplanetary probes and planets themselves. Indeed, since the 
Yukawa-like corrections to the Newtonian accelerations felt by the test particles are usually 
not included in the dynamical force models fit to the observations, an "effective" value of the 
gravitational parameter is, actually, measured, i.e. it is 



Mr 



MoO 



1 + 0C ( 1 + ^ ] 



(12) 



This implies that 



Mmeas 



/'i 



[loo, r > X, 

fioo [1 + 2aexp(-l) 



X, 



(13) 



Mmcas ~ Moo (1+0C), r«l. 
Since a is, of course, expected to be quite small, it is reasonable and adequate to assume 



/'c 



Mr 



(14) 



also for r < X in practical calculations of the perturbative effects of eq. ([9]) (see eq. (|16p below); 
our further analysis will show a-posteriori that this is just the case, given the upper bounds on 
a which will be inferred. Strictly speaking, the use of the measured values /i mca s in those places 
in the formulas in which fj,^ appears would be justified only if X was much smaller than r: this 
would be a fatal restriction because, e.g., //q is routinely measured from interplanetary ranging 
mainly involving the inner planets of the solar system, i.e. one would be forced to consider only 
the case X < 0.38 au = 6 X 10 10 m. 

In view of a first-order perturbative calculation, it is, now, useful to evaluate AUy onto the 
unperturbed Keplerian ellipse and average it over one orbital revolution of the test particle. By 
using the eccentric anomaly E as fast variable, the result is 



(AU- 



v 



exp (-f) 



In 



(15) 



where a, e are the semimajor axis and the eccentricity, respectively, of the orbit of the test 
particle, and Iq(x) is the modified Bessel function of the first kincQ Ik(x) for k = 0. Note that 
eq. (]15j) is an exact result: no approximations have been used either for the orbital configuration 
of the test particle or the size of the scale parameter X. From eq. (|15p it is possible to obtain 
perturbatively the secular precessions of both the pericenter u and the mean anomaly Ai by 
using the Lagrange planetary equations [33] . 
Concerning u, we have 



(d) Y ) 



a 



Moo(l-e 2 )exp(-f) ^ /aeN 
a eX V X J ' 



(16) 



where h(x) is the modified Bessel function of the first kinc0 ifc(x) for k = 1. Notice that eq. 
(|16p agrees with the result obtained by [M] with a different approach. More specifically, 



4 See on the WEB http:/ /mathworld. wolfram.com/ModifiedBesselFunctionoftheFirstKind. html and references 
therein. 

5 See on the WEB http:/ /mathworld. wolfram.com/ModifiedBesselFunctionoftheFirstKind. html and references 
therein. 
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worked out the Yukawa-like pericenter advance per orbit: it can straightforwardly be obtained 
from eq. (fl~6|) by taking the product of (tby) times the orbital period Pb = 2tt/ti = 27r^/a 3 //x 00 . 
The precession of eq. (|16p loses its meaning for e — > since it yields 0/0. Other derivations 
of either the Yukawa-type secular precession of the pericenter or its advance per orbit can be 

found in, e.g., [23 ESI E3 ESI E3 S3 El S3 S3 El S3 S3 S3 SSI S3 ESI ED 153 E3 El ■ Ail of 

them make use of different level of approximations in either the magnitude of the length scale 
X or the orbital configuration of the test particle. From eq. (|16p it is possible to infer 

where 5(Aw) can be thought as the uncertainty in some observationally determined correction 
Am to the standard secular precession of the pericenter for some astronomical system. We will 
use eq. (I17p in Section f3.il and Section T3.2I for the geodetic satellite LAGEOS II and the Moon, 
respectively. 



3 Constraints from laser-ranging 



Let us start to consider the range 10 6 — 10 m corresponding to approximately 1 — 10 Earth's 
radii R c . The traditional constraints on a for X lying in the aforementioned range are depicted 
in Figure 1 of [55] or Figure 4 of [3S], based on Figure 2.13 of [36] and adapted by [56] in his 
Figure 4, Figure 1 of [35], and Figure 1 of [57]: they are of the order of 10 -5 — 5 x 10 -8 . The 
technique with which they have been obtained is described in detail in |36| . It is based on the 
determination of the Earth's gravitational parameter fi c from laser-ranging measurements at 
the altitudes of the LAGEOS satellite and the Moon, and on ground-based measurements of 
the terrestrial gravitational acceleration. 

More specifically, the LAGEOS-Moon constraint of the order of ~ 5 x 10~ 8 is obtained in 
the following way. First, the ratio 



rjA(r h ) - r^Ajru) 
rlA(r h ) + A(ru) 



wher^l 



u) 



L/M 



1 + 0C 



r_L/M\ 



1 + -j-J exp 



%/M\ 

X ) 



[18) 



(19) 



denotes the Newtonian-I- Yukawa acceleration, to be evaluated at distances tl/Mj is theoretically 
computed: eq. (fT8|) is, by construction, independent of //qo. Then, it is compared to the 



analogous ratio 7]'^ computed for the empirically determined values fimcas^ and /iS' of the 
Earth's gravitational parameter at tl and ?"m from Satellite Laser Ranging (SLR) and Lunar 
Laser Ranging (LLR) measurements divided by the square of the Earth's radius R c : in the 
computation of 7? meas we assume Newtonian values for the accelerations ^(^l/m)- Basically, r\ 
is the difference between the values of /i e evaluated at two different distances normalized to the 
average of such two values: in the Newtonian dynamics it vanishes, while in the framework of the 
Yukawa-like deviations from the Newtonian picture 7] is different from zero. The comparison 
between such two determinations of r/ yields upper bounds on a for various values of X. 



3 Here r L / M is a shorthand for designing alternatively either tl or tm- 
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Instead, the LAGEOS-Earth constraint of about 10 5 comes from a comparison between 
the empirical ratio 

• A te „(R c ) - A L {Rc) fon , 

m = MiQ ' (20) 

in which, as usual, standard Newtonian physics is assumed, and the computed one by includ- 
ing the Yukawa term. Notice that A teir (R e ) denotes the acceleration of gravity measured on 
the Earth's surface with ground-based techniques, while A^(R C ) is a shorthand for the ratio 
of 

/iiieas\ i- e - the Earth's gravitational parameter empirically determined with the SLR ob- 
servations to LAGEOS, to the square of the Earth's radius. When r/y is computed, A^(R C ) 
is replaced by the product of the Newtonian-I- Yukawa acceleration ^n+y^l) evaluated at dis- 
tance t*l times (rh/Re) 2 , while At crr (R c ) is replaced by the Newtonian-I- Yukawa acceleration 
evaluated on the Earth's surface A^ + y(R c ). Also r]y is independent, by construction, of fioo. 

Instead, [36] insert /imcas'V-^c i n the denominator of r]y instead of posing ^n+yC^l)!^/^) 2 : 
this is not consistent with all the previously followed line of reasoning, also because in such 
a way /ioo would not be cancelled in rjy. Basically, r/ consists of the difference between the 
accelerations of gravity at the same distance, i.e. on the Earth's surface, measured with dif- 
ferent techniques normalized to the value at the same distance obtained with one of such two 
techniques; again, in Newtonian physics it is expected to be zero, contrary to the Yukawa-type 
case. Note that in rj the value obtained from SLR is extrapolated to the Earth's surface by 
means of the multiplicative scaling factor (r^/R c ) 2 . 

Notice that both for r\ and r) the scale length X has been kept fixed. 



3.1 The perigee of LAGEOS II 

Actually, tighter constraints can be obtained by using the perigee 00 LAGEOS II in connection 
with eq. (TJ6D and eq. flU]). The semi-major axis of LAGEOS II is a = 12, 163 km = 1.2163 x 
10 m = 1.9 R e , while its eccentricity is e = 0.014. 

It is well known that the perigees of the LAGEOS-type satellites are particularly sensitive to 
a host of non-gravitational perturbations which are, thus, a major limiting factor in constraining 
a. As pointed out by Ries et al. [59] in the context of the re-analysis of the earlier tests of the 
general relativistic Lense-Thirring effect with the LAGEOS and LAGEOS II satellites |60| . the 
realistic accuracy in determining the secular perigee precession of LAGEOS II is larger than 
the gravitomagnetic effect itself mainly because of the non-gravitational perturbations affecting 
such an orbital element. From Table 9.7 of |61j it is possible to infer an uncertainty of about 
125 milliarcseconds per year (mas yr _1 in the following). 

Another important source of systematic uncertainty is the mismodeling in the multipoles 
of the expansion of the non-spherical part of the gravitational potential. In particular, the 
fully normalized even (£ = 2, 4, . . .) zonal (?n = 0) Stokes coefficients cause secular perigee 
precessions which acts as superimposed bias. Assessing correctly their mismodeling is crucial 
for a realistic evaluation of the overall uncertainty affecting a. 

There are nowadays several institutions worldwide which almost routinely produce global 
gravity field solutions by analyzing increasing data sets from the dedicated spacecraft-based 
missions like CHAMpH, GRACE, GOCE, along with observations collected from many geodetic 
satellites of LAGEOS-type tracked by the International Laser Ranging Service |62| . Among 

7 For preliminary investigations on the possibility of using the perigee of LAGEOS II to constrain Yukawa-like 
deviations from the Newtonian inverse square law of gravity, see [381 140] . For earlier investigations on the effect 
of a Yukawa-like fifth force on the perigee of LAGEOS-type satellites, see [58] . 

8 It re-entered Earth's atmosphere on September 20, 2010. 
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their products, collected by the International Centre for Global Earth Models (ICGEM) [63], 
the coefficients Ci >m , S^ m of the geopotential stand out. 

For our purposes, it is fundamental to realistically evaluate the uncertainties in the zonals, 
with particular regard to C^o since it induces the largest aliasing secular perigee precession. 
To this aim, we will follow the approach recently described by Wagner and McAdoo [64]. They 
felt the need of comparing not only solutions releasing the mere statistical, formal errors for the 
geopotential coefficients, but also older models yielding calibrated errors, by explicitly requiring 
that the benchmark model must be formally far more accurate than the one to be tested. 
Furthermore, Wagner and McAdoo [64] pointed out that their method could well be applied 
even to solutions not displaying formal errors at all. Wagner and McAdoo |64| remarked that 
the calibration of the errors in a given test model should be made by using reference solutions 
obtained independently: for example, a GRACE-based solution should be compared with, say, 
a GOCE-based solution, as done by Wagner and MacAdoo themselves [64] . Even in such a 
case, care should be taken to avoid that the reference model adopted was not used as a-priori 
background model in producing the models to be tested [64] . 

Here we apply the method by Wagner and McAdoo [64] by choosing the CHAMP-based 
solution AIUB-CHAMP03S [65], the GRACE/GOCE/CHAMP/SLR-based solution GOC02S 
|66j and the GOCE-only solution GOCONS [67] as test models, while we take the wholly inde- 
pendent GRACE/GOCE/LAGEOS-based solution EIGEN-6 [68] as formally superior reference 
model. EIGEN-6 was released by the GeoForschungsZentrum (GFZ), AIUB-CHAMP03S was 
produced by the Astronomisches Institut, Universitat Bern (AIUB), while GOC02S and GO- 
CONS were released by the Gravity Observation Combination (GOCO) consortium. The back- 
ground gravity models adopted for AIUB-CHAMP03S were the pre-CHAMP/GRACE/GOCE 
JGM3 [69] and EGM96 [70] SLR-only solutions, while EIGEN-6 used 6.5 years of LAGEOS 
data from the time span 1 Jan 2003 till 30 June 2009. GOC02S used the GRACE-only solution 
ITG-Grace2010s [71] as background model. GOCONS did not adopt any background model, 
and its data from GOCE cover the period 1 November 2009-17 April 2011. See the discussion in 
|64j about the risk of precluding an unbiased calibration employing an external standard model. 
As requested by Wagner and McAdoo [64] , the formal a-^ 2 of EIGEN-6 is smaller than those 
of the models to be tested by about 1 order of magnitude. Instead of using a single coefficient 

Table 1: First row: uncertainty <5C2,o in the even zonal of degree £ = 2 evaluated according 
to the method by Wagner and McAdoo [04]. EIGEN-6 [68] was assumed as formally superior 
reference model (a^ 2Q = 2 x 10" 13 ), while the solutions AIUB-CHAMP03S [65], GOCONS 
[67] . and GOC02S [66] were taken as test models. The tide system of all the models considered 
is tide-free. Second row: formal a~r< of the aforementioned models, and calibrated a~r< of 

C-2,0 ^2,0 

GIF48, whose tide system is zero-tide. 





AIUB-CHAMP03S (test) 


GOCONS (test) 


GOC02S (test) 


GIF48 


sc 2}0 


3.3 x HT 11 


1.09 x 10~ 10 


1.08 x 1(T 10 




C-2,0 


1.1 x 1CT 11 


1.9 x lcr 11 


4 x 10~ 13 


7.0 x 10" 11 



error scaling factor for the formal sigmas attest of each of the test models, Wagner and McAdoo 
|64] also proposed to average the individual error factors over all the 2£ + l coefficients of degree 
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4xl0" 9 
3.5 x 10" 9 



3 x 1(T 9 
2.5 x 10" 9 

3 

J 2xl0" 9 
o 

1.5 xlO" 9 
1 x 10" 9 
5xlO" 10 

2 4 6 8 10 

A (R e ) 

Figure 1: Upper bound on |oc| for the Earth-LAGEOS II system from the uncertainty in the 
precession of the perigee u) of LAGEOS II as a function of X for 1 R e < X < 10 R e . As far as the 
mismodeling (5(72,0 m the first even zonal harmonic of the Earth's geopotential is concerned, 
we adopted the figures of Table Q] for it, while the uncertainty due to the non- gravitational 
perturbations was taken as large as 125 mas yr _1 |61j . 



I. According to eq. (A13) of [53], one has 



/test/ — 

where H^ m denotes both Ce tTn and Sit >m in the sense that the sum in eq. (I2ip includes all the 
geopotential coefficients of both kinds for a given degree I. Our results are displayed in Table 
[H It can be noticed that the resulting uncertainties 5Ce t o are up to 1 — 3 orders of magnitude 
larger than the formal errors. Incidentally, they are rather close to the calibrated error released 
in the GRACE-only GIF48 model by CSR (5C 2 ,o = 7.0 x 10" 11 ). 

We will adopt our results in Table [1] to assess the zonals-induced bias on the Yukawa- 
induced perigee orbital precession of LAGEOS II. Thus, by using such a figure in eq. (fT71) . it 
is possible to plot the upper bound on a as a function of X: see Figure [TJ It can be noted 
that there is an improvement of 2 — 5 orders of magnitude with respect to Figure 1 of |55j . 
Figure 1 of [35] i and Figure 4 of [32]; see also Figure 1 of [57] . Lucchesi and Peron [72], using 
the approximate analytical results by [40] and a data record 13 yr long for LAGEOS II, claim 
| oc | < 9.9 x 10 -12 at X = 1 R c . However, such a figure should be considered just as a preliminary 
result coming from the statistical errors of the linear fitting of the post-fit numerically integrated 
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A (km) 

1.3 10 4 2.5 10 4 3.8 10 4 5.1 10 4 6.4 10 4 




1 \ 21+1 

' 771=0 



H 



■test 



T7 rcf 



<7- 



ref 

Hi, 



^.test 
77. „ 



(21) 
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2 4 6 8 10 

A (R e ) 



Figure 2: Lower limits on M*/Mp from eq. flHJ for £<0.5 x 10 5 [HE] and the bounds on a of 
Figure CD 



residuals of the perigee since, as pointed out by Lucchesi and Peron themselves [72], it does not 
include an analysis of the systematic errors due to the uncertainties in the zonals' geopotential. 
Lucchesi and Peron [7J used the GRACE-only model EIGEN-GRACE02S [73], which released 
a calibrated a-Q = 5.304 x 10~ . Moreover, Lucchesi and Peron [72] did not explicitly include 
a Yukawa-like fifth force in the dynamical models fitted to the data of LAGEOS II, by assuming 
that it would entirely be accounted for by their numerically integrated residuals of the perigee. 
A similar approach with the LAGEOS satellite was followed by March et al. [74] in constraining 
other non-Newtonian putative effects. 

The bounds on a of Figure [T] yield corresponding constraints on M*/Mp according to eq. 
©. They are depicted in Figured For X = 1 - 1.5 R , it is M*/M P >1.4 - 1.6 x 10 10 , while for 
X = 10 R c we have M*/M P >4-5 x 10 9 . Such results are in striking contrast with M* ~ 10 6 Mp 
inferred by Dvali and Vikman [17] for X = 1 R c . 

3.2 The perigee of the Moon 

Let us, now, consider the range 10 — 100 R e corresponding to 6.4 x 10 7 — 6.4 x 10 8 m; the lunar 
semi-major axis is om = 60 R e = 3.8 x 10 8 m. To this aim, the motion of the Moon is best 
suited to yield tight constraints. 

Reasoning in term of the lunar perigee, the uncertainty in determining its secular precession 
is of the order of 5zb « 0.1 mas yr _1 [75j [76l [77J [751 EH] . The resulting bounds on a are depicted 
in Figure [3j It can be noted that |a| < 3 x 10~ n for X = 60R C = 3.8 x 10 8 m, in substantial 
agreement with Figure 1 of [57] . Figure 1 of [55] and Figure 4 of [39]. Cfr. also with the upper 
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1.2xl0" lu - 



6xl0" n - 



4xlO -11 




20 40 60 80 100 

A (R e ) 

Figure 3: Upper bound on a for the Earth-Moon system from the uncertainty in determining 
the precession of the perigee of the Moon, assumed as large as 0.1 mas yr _1 [T5 j 176 1 \77 \ IT8 | 179] . 
as a function of X for 10 R c < X < 100 R e . 
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20 40 60 80 100 

A (R e ) 

Figure 4: Lower bounds on M*/Mp from eq. (JSj) for £<0.5 x 10 -5 [I] and the lunar bounds on 
a of Figure El 

bound of |oc| < 5 x KT 11 (X = 60 .Re) b>@ Miiller et al. [73 [79]. However, the constraints of 
Figure [3] are orders of magnitude better than those reported in J57J [551 ES] for X 7^ 60i? c . We 
also note that0 [80] obtained |a| < 10" 8 for 1.8 R c < X < 60 i? c . 

In Figure [H we display the corresponding bounds on M*/Mp. While for X = 10 R c it is 
M*/M P ~ 2 x 10 10 , we have M*/M P >3.5 x 10 10 at X ~ 60 i? c , with M*/M P >4.5 x 10 10 at 
X ~ 30 i? e - 

Finally, we remark that using data from the Lunar Reconnaissance Orbiter (LRO) should 
allow to further improve the bounds on |oc| at X = 60 R e |81] , 

4 Discussions and conclusions 

We analytically worked out some orbital effects induced by a hypothetical gravitational fifth 
force of Yukawa-type on the orbital motion of a test particle moving around a central body 
acting as source of the modified gravitational field. We restricted ourselves neither to any 
specific size for the scale length parameter X nor to small values of the orbital eccentricity e of 
the particle. We obtained secular precessions for the pericenter oj. Our results imply the use of 
the modified Bessel functions of the first kind Ik(%), k = 0, 1. 

We used the Satellite Laser Ranging data of the artificial satellite LAGEOS II to put 

9 Miiller et al. [781 179) actually included a Yukawa-type fifth force in the mathematical force models with 
which they analyzed the LLR data. 

10 I thank W.-T. Ni for having pointed out it to me. 
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constraints on the coupling strength parameter a of the putative fifth force for the Earth- 
LAGEOS II system in the range 1 R c < X < 10 R c obtaining 1 x 1CT 10 < |a| max < 4 x 10" 9 . 
The accuracy in determining the Moon's perigee from the Lunar Laser Ranging technique 
allowed us to infer |cc| max ~ 3 x 10~ n for the Earth-Moon pair in the range 15 R c < X < 60 R c . 
From M*/Mp ~ and by using £ < 0.5 x 10 -5 , coming from the overall uncertainty in the 

neutrinic time-of-flight measurements common to both OPERA and ICARUS collaborations, 
our bounds on a correspond to M*/Mp>4 x 10 9 - 1.6 x 10 10 for 1 R c < X < 10 R e , and to 
M*/Mp>3 x 10 10 — 4.5 x 10 10 for 15 R c < X < 60 R e , with the maximum value occurring at 
X ~ 30 R c . Such bounds on M* disagree with M* ~ 10 _6 Mp by Dvali and Vikman. 

Our results are neither necessarily limited to the OPERA neutrino scenario nor to the 
Dvali- Vikman model, being valid for any theoretical scenario yielding an effective long-range 
Yukawa-type correction to the Newtonian inverse-square law. 

As a final, technical remark, we notice that a complementary treatment would require that a 
Yukawa-type extra-force should be explicitly modeled in the softwares used to process the data 
records, and a dedicated solve-or parameter should be estimated in dedicated data reduction 
procedures in which the ad-hoc modified dynamical models are fitted to the existing data sets. 
Such an approach has not yet been implemented in some of the analyses previously described 
like the LAGEOS II one: it is a drawback common to the LAGEOS-based Lense-Thirring 
analyses implemented so far. 
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